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Abstract 

The Schwinger proper-time method is an effective calculation method, ex- 
plicitly gauge invariant and nonperturbative. We make use of this method to 
investigate the radiatively induced Lorentz and CPT- violating effects in quan- 
tum electrodynamics when an axial vector interaction term is introduced in 
the fermionic sector. The induced Lorentz and CPT-violating Chern-Simons 
term coincides with the one obtained using a covariant derivative expansion 
but differs from the result usually obtained in other regularization schemes. 
A possible ambiguity in the approach is also discussed. 
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The question of inducing a Lorentz- and CPT-violating Chern-Simons (CS) term through 
radiative corrections by adding a term of the form tpjs^tp (b^ is a prescribed constant 4- 
vector) to the conventional QED Lagrangian, has been recently addressed by several authors 
|H|-^|. A theoretical investigation on the physical effects of Lorentz and CPT violation 
in electrodynamics has been made by the inclusion of a CS term ^ p k^F v xA p to the 
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Lagrangian density [|],|2|]. The constant 4- vector k^ selects a preferred direction in space- 
time, thus violating both Lorentz invariance and the discrete CPT symmetry. Since up 
to now the observation of distant galaxies does not give any experimental evidence for the 
occurrence of such an effect, the vector k^ is severely constrained and should effectively 



vanish [|IOpT|| . On the other hand, one may ask the question whether such a term can be 
induced through radiative corrections when Lorentz and CPT symmetries are violated in 
other sectors of the theory, for instance in the fermionic one. In this case the stringent 
bounds on k^ would translate into constraints on these sectors. 

The quantum electrodynamics with a Lorentz- and CPT- violating axial vector interac- 
tion term in the fermionic sector is described by the Lagrangian density 

C — ip {ip — eA — m — 75^) ip . (1) 

The full quantum effective action is obtained by integrating out the fermionic fields, 

Z[A] = e iTcft = J V^Vipe l f d4xC = det (ip - eA - m - 75^) , 
T eS = -iTr\n(ip-eA-m-^p) , (2) 

where the trace is taken over both space-time coordinates and spinor indices. 

Several calculations have been carried out to determine the induced CS term, which is 
proportional to b^. While nonperturbative treatments in in the vacuum polarization ten- 
sor and the covariant derivative expansion || yield definite and nonzero values for the 
CS coefficient, the perturbative treatment gives an ambiguous result, which depends on the 
choice of the regularization scheme ||||. Indeed, if one chooses Pauli-Villars regularization 
to enforce gauge invariance for all axial momenta, then the induced mass-independent CS co- 
efficient will vanish JJ. If instead dimensional regularization is used, then different answers 
can be obtained depending on how the matrix 75 is generalized to arbitrary dimensions. 
A universal ambiguous result has also been given in terms of differential regularization J?|. 
Even in terms of the nonperturbative treatment in b^, the calculations based on the Feyn- 
man diagram expansion @J5||| and the non-Feynman diagram method (e.g. the covariant 
derivative expansion) yield different results. 

The purpose of this work is to calculate the Lorentz and CPT- violating Chern-Simons 
term radiatively induced by the axial coupling of the vector b^ in the effective action (0). 
We shall employ the well-known Schwinger proper-time method [Q, which is not only a 
gauge invariant procedure, but also allows us to obtain an explicit solution of the equations 
of motion for the special case of gauge fields which produce a constant field strength tensor 
F^ v . Recently, this approach was employed in combination with the analytic regulariza- 
tion method to investigate the radiatively induced Lorentz and CPT breaking in a constant 



external electromagnetic field |13fl . It was shown that the induced CS coefficient exhibits 
a logarithmic contribution, but this result explicitly contradicts the exact calculation per- 
formed in Ref. |§. As we shall see below, if one follows the original route proposed by 



Schwinger fll2 |, such a nonperturbative and gauge- invariant formulation leads to an induced 
CS current with a finite and definite nonzero value, which is identical to the result obtained 
in the covariant derivative expansion [[J but differs from the result usually obtained in the 
calculation based on the Feynman diagram technique (IHH] . 
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Although the evaluation of the effective action (|2|) is not unique and crucially depends 
on the regularization scheme one adopts, it seems to us that the discrepancies in the re- 
sults obtained for the CS term are not fully due to the chosen regularization. To make our 
point clear, let us first analyze the properties of a regularization dependent ambiguity. As 
is well known, a regularization is a temporary modification of the original theory and is a 
necessary procedure to extract out the divergences so that renormalization can be imple- 
mented. Different regularization schemes actually provide different methods to evaluate a 
quantum correction. For the finite terms in the quantum corrections, the role of a gauge 
invariant regularization is no longer regulating the divergent quantities but rather imposing 
gauge invariance ||. One typical example is provided by the chiral anomaly, which arises 
in any vector gauge invariant regularization scheme such as Pauli-Villars regularization and 
dimensional regularization, but vanishes in any non-invariant regularization method. The 
reason is that these regularization methods preserve the vector gauge symmetry but, simul- 
taneously, violate the axial vector gauge symmetry of the classical theory. On the other 
hand, differential regularization clearly shows that the anomaly manifests itself in both the 
vector and axial vector gauge symmetries. This regularization method does not impose any 
symmetry a priori and to preserve any desired symmetry in the theory one then chooses the 
(undetermined) renormalization scale at the final stage of the calculation [H,7fl. However, 



all the regularization schemes which preserve the same symmetry should give the same fi- 
nite result. If different results are produced by such regularization methods, there must be 
certain inconsistency in the calculation procedure. 

Moreover, it has recently been claimed that if one uses a covariant derivative expansion 
to calculate the coefficient of the radiatively induced CS term, then there is a discrepancy 
with the results previously obtained ||,]5]|§ by using the widely accepted approach of the 
chiral anomaly || . It was explicitly shown in Ref . || that when compared with the standard 
perturbative expansion in powers of the gauge field A^, the one- loop vacuum polarization 
tensor in the covariant derivative expansion has an additional term, whose surface term 
gives a non-vanishing contribution to the usual result and which is precisely the source of 
the discrepancy. Since the essence of the covariant derivative expansion is to develop the 
local effective Lagrangian in a series of powers of ILj = id^ — eA^, rather than id^ and A^ 
separately [|l5|,[n|, this formalism is explicitly gauge covariant. In this sense the covariant 



derivative expansion and the Schwinger proper-time method share the same features, and 
this is another motivation for us to employ the latter method in our calculation. 
Let us first write the trace in Eq. (fj) in the following equivalent form, 

Tr In (ift — eA — m — 75$) = Tr In (id — eA — m) — [ dzTi \ — -75$ | . (3) 

jo \i<p — efl — m — zj 5 p J 

Then the quantum effective action can be rewritten as 

r* = r2 + r2, (4) 



where 



r (o) 

1 cff 



r% = i dzTr - 75 . (5) 
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The radiatively induced CS term comes only from the T^g part. Moreover, since Tqs oc b 
we can neglect the dependence of b in the denominator of the trace in T^. Thus we have, 



cs 



iTr 



io — €m — m 



iTr 



II — m 



75 1 



(6) 



with the commutation relations [x^, IIj,] = —ig^ v > [n M , U u ] = —ieF^. 

To calculate the radiatively induced CS term we make use of the Schwinger proper-time 



(7) 



method ||12|| . We begin by writing Eq. (|6|) as 

Tcs = i J d 4 xTr (^( x \ j^ ^ m \ x )^^j = i J d^xTr [G(x,x')j 5: 

where the function G(x, x') is defined by 

(g-m) G(x, x') = 5 (4) (a; - x') , (8) 
and the limit x' — > x is performed by taking the average of the terms obtained by letting 



x' — > x from the future and from the past |T!|. In operator form, 

1 „ , 1 1 



G(x,x') = (x\G\x'), G 



II — m 1/t 2 — m 2 Tjl 2 — m? 

Therefore we have two equivalent integral representations for the operator G 

dse im2s U{s) fll + m), 



(9) 



G = -i $ + m) / dse tms U(s) = - 

where U(s) = e~ lHs and 

U=^ = ll 2 - e -o-F, c-F = o ilv F»\ a ll v= l -[l tl ,lv}- 
Eqs. (|^) and (|TDD show that the calculation of G(x,x') lies in the evaluation of 

(x'\U(s)\x") = (x'(s)\x"(0)) . 



(10) 



(11) 



(12) 



The operator U(s) can be regarded as an evolution operator of a system governed by the 
"Hamiltonian" 7i in the "time" s. The Heisenberg equations of motion for the operators 



x M (s) and TL^s) read 



.dx. 



ds 

ft. 

ds 



-2iU 



[n M , H] = -ie (F^IT + Tl v F^y) - -iea up d»F up 
-lieF^lY + ed u F uu - Uea up d u F 



H± up- 



pu g. 

Correspondingly, the transformation function satisfies the differential equations 



(13) 
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i^(x'(s)\x"(Q)) = (x'(s)\H\x"(Q)), 



d 

-i-z eAJx") 



(x'(s)\x"(0)) = (x\s)\U fl (s)\x"(0)), 
(x'(s)\x"(0)) = (x'( S )|n,(0)|x"(0)), 



with the initial condition 



(x'(s)\x"(0))\ s ^ = 6W(x'-x" 



(14) 



(15) 



We shall now solve the Heisenberg equations (|T3| ) for the special case of a constant field 
strength tensor, i.e. we choose to be constant, = —F fJiU x u /2. This will allow us to 
calculate the transformation function (x' (s)\x" (0)) in an exact analytical form. In this case 
the proper-time dynamical equations (|T3|) are simplified to 



ds 



-2n 



ds 



or in matrix notation, 



ds ' ds 



-2eF^rT 



-2eFIT . 



The formal solutions can be easily written out, 



(16) 



(17) 



U(s) 



-2eFs 



n(o) , x(s)-x(o) 



eF 



e -2eFs _ j 



n(o) 



To get the transformation function, we express n(s) in terms of x(s), 

11(0) = --eFe^sinlr^eFs) [x(s) - x{0)} , 
n(s) = — eFe-^'sinh-^eFs) [rr(s) - x(0)] 
= - [z(s) - x(0)} -eFe eFs sinlr 1 (eFs) . 

Therefore, 



(18) 



(19) 



n 2 (s) = [x(s) - x(0)] X [x(s) - z(0)] 
K = -e 2 F 2 sinh" 2 (eFs). 



(20) 



Using the commutation relation 



x(s),x{0)} = [x(0) + (eF)" 1 (e~ 2eFs - l) 11(0), x(0) 



KeF)- 1 ( 



e 2eFs - 1 



(21) 



we can write 
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n 2 (s) = x(s)Kx(s) - 2x(s)Kx(0) + x(0)Kx(0) - -Tr [eFcoth (eFs) 



(22) 



where the trace extends over space-time indices. 

Taking into account Eq. fl22|) , the first differential equation in Eq. ( |i~4|) becomes 

i-^-(x'(s)\x"(0)) = {--a • F + {x' — x")K(x' - x") 
as 12 



which has the solution 



where 



(x'{s)\x"(0)) = C{x',x")e 



L(s) = -Trln 
v ; 2 



--Tr [eF coth (eFs)] \ (x (s)\x" (0)) , 



I ^'\ P -L(s) s -2 e {{x' -x")eF coth(eFs){x' -x") ^ea-Fs 



(eFs) 1 sinh(eFs) 



(23) 



(24) 



(25) 



The coefficient function C(x',x") can be determined from the second and third differential 
equations in (0). Since from Eq. QT§D we have, 



(x'(s) |n(0) \x"(0)) = - - [eF coth (eFs) + eF) (x' - x") (x'(s) \x"(0)) 
(x'(s)\Tl(s) \x"(0)) = -- [eF coth (eFs) - eF] (x' - x") (x'(s) \x"(0)} 



(26) 



then Eqs. fll4|), fl24]) and (^) lead to the following differential equations for C(x', x") 

C(x', x") = , 



i^-eA,(x')- e -F, v (x'-x"Y 



-^-eA»(x») + e -F» w (x'-xr 
The solution of these equations is given by 

C(x', x") = C$(x\ x") , $(x\ x") = exp 



C(x', x' 1 







(27) 



ie / dx^A^x) 



(28) 



where we have chosen the integration path to be a straight line connecting x' and x", since 
the integral in the exponent of Eq. (|28D is independent of the path. The overall constant C 
can be fixed by the initial condition ( |I~5"D and we find 



(29) 



C 



(4tt) 2 



Thus we finally get the transformation function 



(x'(s)\x"(0)) 



(4tt) : 



$(x', x ") e - L ^) s -^ e i{x'-x")eFcoth{eFs)(x'-x") e iea-Fs ^ 



(30) 
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and the corresponding Green function 



G(x',x") = -i \ dse m2s [Y(x'(s)\U ll (s)\x"(0)} +m(x'(s)\x"(0))] 
= -if° dse im2s [(x'(s)\U fl (0)\x , '(0)}^ + m(x' (s)\x" (0))] . 



(31) 



We are now in position to obtain the radiatively induced CS term (0). Using the trace 



relation Tr 



0, we have 



Tr [G(x,x'h0] = -tTr T6 0f / ds e im " s (x(s)\U,(s)\x'(0)) 



Tr 



— iTr 



il^ + w^bv) / ^e jm2s (x( S )|n M ( S )|x'(0)) 



7^/ dse ims (x(s)\U^0)\x'(0)) 



Tr 



(ij^ + jsa^K) / dse im2s {x(s)\UJ0)\x'(0)) 



(32) 



On averaging these two equivalent expressions we find that 



Tr[G(x,x')j 5 #\ =-OTr| 75 ds e im2s {x(s)\ l - [U,(s) - II M (0)] \x'(0)) 
+ KTr Ua^ [° dse m2s (x(s)\l [Tl^s) + iyo)] |x'(0)> 



(33) 



According to Eqs. P 



P„(«) - n„(o)] |z'(0)> = |f m „(x - xO^x^lx^O)) , 
<x(s)|± [n M ( S ) + n M (0)] \x'(0)) = —F^ [coth(eF S )] Q/3 (x - x%(x(s)\x'(0)) . 



(34) 



Thus we obtain 



Tr [G(x, x')^ 5 j§\ = --M^Tr 



^ v I dse im ' s [coth{eFs)] al3 (x - 3/) p (x(s)\a/(0)) 



re 



-b„F^(x - x') w Tr 



75/ dse im s (x(s)\x'(0)) 



(35) 



Using Eq. ( |3~0| ) we can write then 



Tr[G(x } x'h^ 



:§(x,x')b y 



ds 



327T 2 

x {i^Tr 
- iF» v {x - x%Tr 



gim 2 Sg— L(s) j(x— x')eF coth(eFs)(x— x') 



oo S 



75e 



coth(eFs)] a/3 (x - x') p 



(36) 



It turns out to be more convenient to calculate the ground-state current in the presence 
of the background field A^x), 



W*)) 



^r C s 



(37) 



From Eqs. (|7|) and (|36|) , and after deforming the integration path by the substitution s — > —is 
we get the nonperturbative expression for the CS current: 



X-X')^(x,x')b v f° ^ e m\-L(s) e -\(x-x >)eF cot(eFs)( X - X >) 



327T 2 

x {-F pa Tr 



[cot(eFs)] a/3 (ar-x')j8 

#<r-Fs" 



75 e 



(38) 



where the operator L(s) is given by Eq. (p5|) after the corresponding substitution s — ► — zs. 

To evaluate the CS current, the following operations should be performed. First we 
notice that since the gauge-invariant quantity e~ L ^ is a determinant, it can be evaluated 
from the eigenvalues of the matrix F, which can be easily found with the assistance of the 
relations: 



pwp 



F w F pv - F iip F pu = 25; T, 



(39) 
(40) 



where F^ v = \e^ xp F Xp is the dual field strength tensor, T = \F^ V F^ U = -\F^F pu 



4 v 



E 2 — B 2 ) and Q = -F^F^ = — E • B. We iterate the eigenvalue equations: 



F^ v = Xipp , F^tp" = -\pi>ti > 



to obtain 



/^./•'"v,, a 2 ^, f^f^ = ^g 2 ^ . 



The identity ( f40|) yields then the equation 

A 4 + 2FA 2 - g 2 = 
which has four solutions, ±Ai and ±A2, 



Ai 



% 

< 2 = T2 



Vf+igtVF- Jo 



(41) 



(42) 



(43) 



(44) 



The quantity e L ^ can now be expressed in terms of these eigenvalues, 



-L(s) 



2 ^1 



A, 



(es) 2 £ 



sin(eAis) sin(eA 2 s) ImfcoshesX] 



(45) 



where X = ^J2{F + ig). 

Next we need to evaluate the Dirac traces. It is easy to show that the following matrix 
decomposition holds: 
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\cr-Fs 



where the coefficients q are given by 

ci = Recosh(esX) , c 3 = Imcosh(esX) , 

c 2 = Re [sinh(esX)/2X] , c 4 = Im [sinh(esX)/2X] . 

In the expression ([38]) for the CS current the following Dirac traces appear: 

Tr{ l5 e^ Fs } = 4zc 3 , 
Tr{7 5 (7^ef <7 - Fs } = 8i [c 2 F pv + c 4 F^) . 

In deriving the above formulas we have used the relations: 

Tr 75 = Tr o pv = Tr^a^} = , 
Ti{a pu a a(3 } = A(g pa g uP - g pP g ua ) , 
Tr {7 5 <x p V Q/3 } = Ait pual3 . 



(46) 



(47) 

(48) 
(49) 



(50) 



Substituting Eqs. 



and (f49 ) into Eq. (^) we obtain 



( Jp{x)) = (x - x%$(x, x')fy,£ y c?s 
i^pa (c 2 F" v + c 4 F^) [cot(eFs) 



Im[cosh(esX)] 



^(x-x')eF cot(eFs)(x—x') 



I a/3 



C 3 



(x - x% - ^(x - x') p 



(51) 



For our purposes it is sufficient to evaluate Eq. ( |5TD in the weak field approximation. In this 

case 



Im[cosh(esX)] ~ {es) 2 Q , eFcot(eFs) : 
es 

c 2 - — , c 3 ~ , c 4 ~ , 
and therefore the leading term in the CS current (|5i~|) reads as 

/ m 2 K 1 (j— m 2 (x — x') 2 ] 1 ^ 2 ) 



1 



(52) 



(■/*.(*)> 



2vr 2 



-$(x, x 



-m 2 (x — x') 2 ] 



(53) 



where if 1(2) is the first-order modified Bessel function defined through its integral represen- 
tation 2 y /zK 1 (2 y /z) = f™ds s~ 2 e" 52 - 1 /*. Finally, making use of the following formulae 



^(X, X ) Ijj'^-j; 



Ki ([— m 2 (x — x') 

1 •> 



— m x — x 



21 1/2 



m 2 (x — x') 



A2 



(X X )/i(X X 



x — X 



A2 



d 4 x 



/ (X X )^i(X X )^y 



x — X 



A2 



(54) 
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we obtain 



(Mx)) = —F^bu = ~^- 2 e^b u F Xp . (55) 

By functionally integrating over the gauge field A p (cf. Eq. ([J7|)) we get the radiatively 
induced Lorentz and CPT violating CS action 



e 2 



cs 



16tt 



J d A xe^ x %A„F Xp . (56) 



Eq. (|56|) shows that the Schwinger proper-time formulation yields a finite and nonzero 
value for the radiatively induced CS term. We note that this result coincides with the one 
obtained in Ref. || using a covariant derivative expansion, but it differs from the result 
of Refs. jiHH], wnere a coefficient 3/(327r 2 ) is obtained instead of l/(167r 2 ). As discussed 
before, this discrepancy is not due to the ambiguity in the choice of a regularization scheme. 
The covariant derivative expansion shows clearly that the vacuum polarization tensor re- 
ceives a non-Feynman diagram contribution ||, which is the origin of the discrepancy in 
the results for the Chern-Simons coefficient. The Schwinger constant field approximation 
method has a common feature with the covariant derivative expansion, namely, the (axial 
vector gauge) anomalous effective action is constructed from (vector) gauge invariant (or co- 
variant) quantities. This is probably the reason why these two methods lead to an identical 
result. 

We should however emphasize that the results given in Eqs. (|55|) and (|56"D have a potential 



ambiguity. This can be seen from the last limit in Eq. (|54|) , which is, in a strict mathematical 
sense, directional dependent, i.e. it may take distinct values as x p approaches to zero from 
different directions. In general, \im x ^ x fl x u /x 2 = Cg pu with C being an arbitrary constant. 
To ensure the consistency of the trace in four dimensions, we choose C = 1/4. It is in this 
sense that the value of the radiatively induced Lorentz and CPT violation is ambiguous in 
our approach. From a theoretical point of view, we find the calculations performed using 
the covariant derivative expansion and the Schwinger proper-time method more appealing. 
Nevertheless, if the radiatively induced Chern-Simons term has any physical observable 
effect, it is the comparison with experiment which will fix such ambiguities and ultimately 
resolve the discrepancies among various results. 
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